Semileptonic and radiative decays of the Be meson in light-front quark model 
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We investigate the exclusive semileptonic Be {D,ric, B , Bs)£vi, rjb B^lvill — e,jj,,T) decays 
using the light-front quark model constrained by the variational principle for the QCD motivated 
effective Hamiltonian. The form factors f+{q^) and f-{<f) are obtained from the analytic contin- 
uation method in the g+ = frame. While the form factor f+{q^) is free from the zero-mode, the 
form factor f-{q^) is not free from the zero-mode in the q'^ = frame. We quantify the zero-mode 
contributions to f-{q^) for various semileptonic Be decays. Using our effective method to relate 
the non-wave function vertex to the light-front valence wave function, we incorporate the zero-mode 
contribution as a convolution of zero-mode operator with the initial and final state wave functions. 
Our results are then compared to the available experimental data and the results from other the- 
oretical approaches. Since the prediction on the magnetic dipole B* —> Be + y decay turns out to 
be very sensitive to the mass difference between B* and Be mesons, the decay width r(_B* Bej) 
may help in determining the mass of Be experimentally. Furthermore, we compare the results from 
the harmonic oscillator potential and the linear potential and identify the decay processes that are 
sensitive to the choice of confining potential. From the future experimental data on these sensitive 
processes, one may obtain more realistic information on the potential between quark and antiquark 
in the heavy meson system. 



> 

p 

rn 
O 

o 



X 



I. INTRODUCTION 

The exclusive semileptonic decay processes of heavy 
mesons generated a great excitement not only in ex- 
tracting the most accurate values of Cabbibo-Kobayashi- 
Maskawa(CKM) matrix elements but also in testing di- 
verse theoretical approaches to describe the internal 
structure of hadrons. The great virtue of semileptonic de- 
cay processes is that the effects of the strong interaction 
can be separated from the effects of the weak interaction 
into a set of Lorcntz-invariant form factors, i.e., the essen- 
tial informations of the strongly interacting quark/gluon 
structure inside hadrons. Thus, the theoretical problem 
associated with analyzing semileptonic decay processes is 
essentially that of calculating the weak form factors. 

In particular, along with the experimental study 
planned both at the Tevatron and at the Large Hadron 
Collider(LHC), the study of the Be meson has been very 
interesting due to its outstanding feature; i.e., the Be me- 
son is the lowest bound state of two heavy (6, c) quarks 
with different flavors. Because of the fact that the Be 
meson carries the flavor explicitly, not like the symmet- 
ric heavy quarkonium (66, cc) states, there is no giuon 
or photon annihilation via strong interaction or electro- 
magnetic interaction. It can decay only via weak inter- 
action. Since both 6-and c-quarks forming the Be meson 
are heavy, the Be meson can decay appreciably not only 
through the b q{q = c, u) transition with c quark being 
a spectator but also through the c q{q — s, d) tran- 
sition with b quark being a spectator. The former tran- 
sitions correspond to the semileptonic decays to rje and 
D mesons, while the latter transitions correspond to the 
decays to B^ and B mesons. The latter transitions are 
governed typically by much larger CKM matrix element; 
e.g., \Ves\ ~ 1 for Be ^ BsfMeil - e,n), vs. \Veb\ 0.04 
for Be rietvi{£ = e, /i, r). For this reason, although the 



phase space in c ^ s,d transitions is much smaller than 
that in 6 —> c, u transitions, the c-quark decays provide 
about ~ 70% to the decay width of Be- The 6-quark 
decays and weak annihilation add about 20% and 10%, 
respectively [l|. This indicates that both 6-and c-quark 
decay processes contribute on a comparable footing to 
the Be decay width. 

There are many theoretical approaches to the calcula- 
tion of exclusive Be semileptonic decay modes. Although 
we may not be able to list them alL we may note here the 
following works: QCD sum rules 0, 0, S 13 , the relativis- 
tic quark model [1^ 0, based on an effective Lagrangian 
describing the coupling of hadrons to their constituent 
quarks, the quasipotential approach to the relativistic 
quark model E&IS: 

the instantaneous nonrelativistic 
approach to the Bethe-Salpeter(BS) equation [ll[, the 
relativistic quark model based on the BS equation [l^, 
[l3| . the QCD relativistic potential model ^14|, the rela- 
tivistic quark- meson model [15*1, the nonrelativistic quark 
model 16], the covariant light-fro nt q uark model [13], 
and the constituent quark model O [S l20l [2ll us- 
ing BSW(Bauer, Stech, and Wirbel) model and 
ISGW(Isgur, Scora, Grinstein, and Wise) model [23]. 

The purpose of this paper is to extend our light-front 
quark model(LFQM) [H, [H [13, [H, based on the 
QCD-motivated effective LF Hamiltonian to calculate 
the hadronic form factors and decay widths for the ex- 
clusive semileptonic Be Pivi{P = D, Dg, B, Bg) and 
77fc — > Beii'e decays and the magnetic dipole Bl Bcj 
transition. In our previous LFQM analysis [2J, [2^, [2y, 
l27i . l2g . [29l . we have analyzed the meson mass spec- 
tra |24l . [25| and various exclusive processes of the ground 
state pseudoscalar(P) and vector(U) mesons such as the 
P ^ P semileptonic heavy/light meson decays |2^ [26|, 
the rare B — > Ki'^i^ decays [13], and the magnetic 
dipole transitions of the low- lying heavy/light pseu- 
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doscalar/ vector mesons [13, [H, [2^. In those analyses, 
we found a good agreement with the experimental data. 
However, since we didn't analyze the Be and B* mesons 
yet, we shall extend our LFQM to predict the masses and 
the decay constants of Be and B* mesons as well as the 
above mentioned exclusive decays of Be and B* mesons. 

Our LFQM [li, H [13, S [H analysis in this work 
has several salient features: (1) We have implemented 
the variational principle to the QCD motivated effective 
LF Hamiltonian to enable us to analyze the meson mass 
spectra and to find optimized model parameters. The 
present investigation further constrains the phenomeno- 
logical parameters and extends the applicability of our 
LFQM to the wider range of hadronic phenomena. (2) 
We have performed the analytical continuation from the 
spacelike region to the physical timelike region to obtain 
the weak form factor f+{q^) for the exclusive semilep- 
tonic decays between the two pseudoscalar mesons as well 
as to obtain the decay form factors Fvp{q^) for V Pj* 
transitions. The Drell-Yan-West((7+ = g'^ + g'' = 0) frame 
(i.e., = — q5_ < 0) is useful because only the valence 
contributions are needed unless the zero-mode contribu- 
tion exists. 

The form factor /+((?^) can be obtained just from the 
valence contribution in the = frame without en- 
countering the zero-mode contribution [30l |. However, 
the form factor f-{q^) receives the higher Fock state 
contribution(i.e., the zero-mode in the q'^ — frame or 
the nonvalence contribution in the g"*" > frame) within 
the framework of LF quantization. Thus, it is necessary 
to include either the zero-mode contribution(if working 
in the g"*" = frame) or the nonvalence contribution(if 
working in the g"*" > frame) to obtain the form fac- 
tor f-{q^)- In this work, we utilize our effective method 
presented in to express the zero-mode contribution 
as a convolution of zero-mode operator that we find in 
this work with the initial and final state LF wave func- 
tions. In this way, we calculate the form factor f-{q^) in 
the q^ — Q frame with the perpendicular components of 
the currents and discuss the LF covariance of f-{q^) in 
the valence region by analyzing the covariant BS model 
and the LF covariant analysis described by Jaus [3l| . We 
also estimate the zero- mode contributions to the f-{q^) 
for various semileptonic Be decays in our LFQM. 

The paper is organized as follows. In Sec. II, we dis- 
cuss the P ^ P semileptonic decays using an exactly 
solvable model based on the covariant BS model of (3-1-1)- 
dimensional fermion field theory. We explicitly show the 
equivalence between the results obtained by the mani- 
festly covariant method and the LF method in the (7+ = 
frame. The extraction of the zero-mode contribution to 
f-^{q^) in the g+ = frame and the effective inclusion 
of the zero-mode in the valence region are discussed. In 
Sec. HI, we briefiy describe the formulation of our LFQM 
and the procedure of fixing the model parameters using 
the variational principle for the QCD motivated effective 
Hamiltonian. The masses and decay constants of the B* 
and Be mesons are predicted and compared with the data 



as well as other theoretical model predictions. The dis- 
tribution amplitudes(DAs) for the heavy-ffavored mesons 
such as D,r]e, B, Bs, Be and iji, are also obtained in this 
section. In Sec. IV, we calculate the weak form factors 
f+{q^) and f~{q^) in the q+ — Q frame using the plus 
and perpendicular components of the currents, respec- 
tively. The zero-mode contribution to the form factor 
J~{q^) is also discussed. In Sec. V, the decay form factor 
FbibAi^) foi' Bl — > Be"f* transition and the decay 
width for B* Be"f are presented. The coupling con- 
stant gB»B^ needed for the calculation of the decay width 
for B* — > -Bc7 is determined in the limit q^ — s- 0, i.e., 
gs^B^ = FB'sS'f = 0)- For the numerical calculation 
of the semileptonic and radiative decays, the form factors 
I±{q^) for the semileptonic decays and Fb-'bAi^) fo'' the 
B* — > -Bc7* transition are analytically continued to the 
timelike g^ > region by changing to —q^ in the 
form factor. In Sec. VI, our numerical results (i.e., the 
form factors and decay rates for Be {D,rie, B, Bs)ii'i, 
rib — > Beii'e, and B* — > Bej^*^ decays) are presented and 
compared with the experimental data as well as other 
theoretical results. Summary and discussion follow in 
Sec.VII. 



II. P SEMILEPTONIC DECAYS IN 

COVARIANT BETHE-SALPETER MODEL 

A. Manifestly covariant calculation 

The amplitude A for a semileptonic decay of a meson 
Qiq with the four-momentum Pi and the mass Afi into 
another meson Q2q with the four- momentum P2 and the 
mass M2 is given by 

A=^Vq^q^L,H'^, (1) 

where Gf is the Fermi constant, Vq_^q^ is the relevant 
CKM mixing matrix element, is the lepton current 

L^^u^,'yf,{l~j'')ve, (2) 
and H^^ is the hadron current 

i7^ = (P2,e|(^^-^^)|Pi). (3) 

Here, e is the polarization of the daughter meson and 
and are the vector and axial vector currents, re- 
spectively. If the final state is pseudoscalar, the hadron 
current can be decomposed as follows: 

{P2\A'^\Pi) = 0, 

{P2\V^\Pl) = .f+{q^){Pi+P2Y + f-{q^)q^, (4) 

where q^ = (Pi — P^)^ is the four-momentum transfer 
to the lepton pair(^t'^) and mj < q"^ < (Mi — A'h)^- 
Sometimes it is useful to express the matrix element of 
the vector current in terms of f+{q^) and fa{q^), which 
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correspond to the transition amplitudes with 1~ and 0"*" 
spin-parity quantum numbers in the center of mass of 
the lepton pair, respectively. They satisfy the following 
relation: 



f-{q' 



(5) 



Including the nonzero lepton mass, the differential de- 
cay rate for the exclusive 0^ O^iiy^ process is given 



x|[AV)P(l + ^jl/+(9')l= 

where K{q^) is the kinematic factor given by 



(6) 



1 



2Mi 



{Ml+Ml-q^Y-AMllvq. (7) 



The solvable model, based on the covariant Bcthe- 
Salpeter(BS) model of (3 -I- l)-dimensional fcrmion field 
theory (ssl . [3^ . [ssj, enables us to derive the transition 
form factors between two pseudoscalar mesons explicitly. 
The matrix element = (P2|V^^|-Pi) in this BS model 
is given by 



(8) 



where gi and g2 are the normalization factors which can 
be fixed by requiring both charge form factors of pseu- 
doscalar mesons to be unity at zero momentum transfer, 
respectively. To regularize the covariant fermion triangle- 
loop in (3-1-1) dimensions, we replace the point gauge- 
boson vertex 7^(1 — 75) by a non-local (smeared) gauge- 
boson vertex (Ai^/TVa J7''(l - 75)(A2^/A'"a2), where 
TVaj = Pi - Ai^ + ie and N\.^ ^ pI ~ ^-2^ + ie, and thus 
the factor (AiA2)^ appears in the normalization factor. 
Ai and A2 play the role of momentum cut-offs similar to 
the Pauli-Villars regularization [H, [s^l . The rest of the 
denominators in Eq. i.e., NiNqN2, are coming from 





= pI 


— rrii^ 


-f ie 








-ie, 


N2 


= pI 


- m2^ - 


-f ie 



the intermediate fermion propagators in the triangle loop 
diagram and are given by 



(9) 



where mi, ruq, and m2 are the masses of the constituents 
carrying the intermediate four-momenta pi = Pi — k, k, 
and P2 — P2 ~ k, respectively. Furthermore, the trace 
term in Eq. ([8]), S'^, is given by 

5^ - Tr[75(^i + mi)7''(:^2 + m2)75(- ^ + m,-)] 
= A^k- P2-k^ + m2mq)Pi' 
+4:{k ■ Pi-k^ + mimq)P^ 
-|-4(fc^ — Pi ■ P2 — mirriq — m2'mq -\- mim2)k^ . 

(10) 

We then decompose the product of five denominators 
given in Eq. ([S]) as follows: 

1 1 



iVAiiVliV,-A^2A^A, (Ai2-TOi2)(A22_TO22) 

1/1 1 \ / 1 1 



Nq\NA, NJ\Na, N2 



(11) 



Once we reduce the five propagators into a sum of 
terms containing three propagators using Eq. pip , we 
use the Feynman parametrization for the three propaga- 
tors, e.g.. 



1 



dx I dy 

NiNqN2 Jo Jo 



[Nq + {Ni - Nq)x + (Ar2 - Nq)yY 



(12) 



We then make a Wick rotation of Eq. ([8]) in _D-dimensions 
to regularize the integral, since otherwise one looses the 
logarithmically divergent terms in Eq. ([H]). Following the 
above procedure, we finally obtain the Lorentz-invariant 
form factors f+{q^) and f-{q^) as follows: 
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Ml') 



f-ii') = 



N 



87r2(Ai^-mi2)(A2^-m22)7o 



dx 



l-x 



dy\ [3{x + y) - 4] In 



CAiA2C'min 



(1 — .T — y) {xM-^ + yAfj) + a;?;(2 — a; — y)^/ + (x + y){mim2 — mimg — m2'mq) + m^(mi + 7712) 



87r2(Ai^-mi2)(A2"-TO22)7o 



dy<^ 3(x - y) In 



CAim2C'miA2 
Cai A2 C^nira-z 



{yMl - xMl) 



+ {x^ - y'^){xMl + yM|) - xy(a; - y)q^ + (x - y){mini2 - toito^ - m2mq) + m^(m2 - mi) 



where iV = 51^2 A|A^ and C = (I/CA1A2 - l/CAim2 - l/CmiA2 + 1/Cmim2) with 



CA1A2 


= (1 


— X — 


y)(xAf2- 


H yM|) - 


h xyg^ 


- (xA^ ^ 




(l-x- 


CAim2 


= (1 


— X — 


y)(xAf2- 


1- yAf2) - 


h xyg^ 


- (xA^ ^ 


-yml) - 


(l-x 




= (1 


— X — 


y)(xMi2- 


K yAf|) - 


h xyq^ 


— {xrrii - 


fyA^)- 


(l-x 


Crnim2 


= (1 


— X — 


y)(xAfi^- 


K yAf|) - 


h xyg^ 


— (XTO^ - 


f yml) - 


- (1 - X 



y)m?. 



9' 
^1- 



c 



(13) 



Note that the logarithmic terms in f+{q^) and f-{q^) are obtained from the dimensional regularization with the 
rotation. 



(14) 
Wick 



B. Light-front calculation 



Performing the LF calculation of Eq. ([5]) in the (7+ = frame in parallel with the manifestly covariant calculation, 
we shall use the plus and perpendicular components of the currents to obtain the form factors f+{q^) and f-{q^), 
respectively. That is, in the q~^ = frame, one obtains the relations between the current matrix elements and the 
weak form factors as follows 



2P^ 



(15) 



The LF calculation for the trace term in Eq. (jTU]) can be separated into the on-shell propagating part 5*^,-^ and the 
instantaneous part S'j'^^^ via 



as 



where 



and 



CM — CM I CM 
— '^on + '-'inst' 



CM 
"-"on 



Plon(P2on ' fcon) - ^onblon ' P2on) + P2on(Plon -^011)+ m2m^p5'on + "llTOgP^^^ + mi 7712 fc^,, 



P2on^<ii - PtonKn + ff^(P2on ' fcon + m2mg) 



+ 2(P2 ~P2on) 



+2{k- - fc- ) 



Plon^mi - PtonKn + 9^^ {Plon ' + mim,) 

PlonP2on + PlonP2o„ " ^''^(Plon -^2011 " mim2) 



+ 2g'"+fc+„(Pl -Pio„)(j32 -P2on)- 



(16) 

(17) 
(18) 



(19) 
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Note that the subscript (on) denotes the on-mass-shell {p = m ) quark momentum, i.e., p = p^^ = (m + p\)/p 
The traces in Eqs. and are then obtained as 



4P+ 

S+ = -^ik^-k'^+AiA2), 



(20) 



for the plus component of the currents and 
-2k I 



x{l — x) 



2k^ ■ k'^ + (1 — x){q~i + ml + 1712) + Ix + 2x{l — x){raimq + m2mq — mim2) 



-(kl+A\), 



x{l — x) 



"-■inst — ^-^ 1 



{Pi - Pio„)k'± + {P2 - P2on)k± + x{k - fcon)(2k± + qj 



(21) 



for the perpendicular components of the currents, where k'j^ = kj^ + (1 — a;)q^ and Ai = (1 — x)nii + xruq (i = 1, 2). 

As one can see from Eqs. (|20p and (|2ip . the perpendicular components of the currents receive instantaneous contri- 
butions while the plus component of the currents does not receive them. Especially, the absence of the instantaneous 
contributions to the plus current indicates that there is no zero-mode contribution to the hadronic matrix element of 
the plus current. 



1. Valence contribution 



In the valence region < A;+ < , the pole k ~ k^^ — {]<.']_ + m? — ie)/k^ (i.e., the spectator quark) is located 
in the lower half of the complex A;^-plane. Thus, the Cauchy integration formula for the k~ integral in Eq. ([8]) gives 

N dx 



N r dr r 



(22) 



where S>^^^ = + S^^Jk- = fc" ) and 



1 - X 



k'l + ml k' 



1 - X 



(23) 



and Ml^ = M^{mi Ai), M^^ ^ Ml^{m2 A2) with k'_L = k_L + (1 - a;)q^. The LF vertex functions xi and X2 
are given by 



Xi(x,k^_) 



1 



-, x2(a;,k'_L) 



1 



From Eqs. ([^0)1 and (PT|) . we obtain the valence contribution to -M^^^i and A^^^; as follows 



(24) 



m: 



47^3 
N 



r (13^ J d^^±Xi{x,k±)x2{x,k'j_){k±-k'^+AiA2), 



dx 



167r3 Jo (l-x) 



X 1 (2; , ) X2 (a; , k' ^ ) 5*4; , 



where 



Svai - -2kj_ 



Aff + Mi + qi - (mi - mqf - (mz -mg)^ + (mi - ma)^ 



2qj 



(1 - x)Ml + xMq - (mi - m,)^ 



(25) 



(26) 
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From Eqs. ((T5)l and p5)) . we get the LF valence contributions to f+{q^) and f-{q^) as follows 

N rlr r 



X — (1 — x)Afi + (to2 — rnq)Ai — mg(mi — TOg) H ^-^^[Aff + M| — 2(toi — mq){m2 — m,)] 



/™_ ™ , /I ™ /™ ™ 1 1 — z,yiiii — inq)(iii2 — iiiq)\ f 

(27) 

We note that the form factors in Eq. (|27p obtained in the spacelike region using the — frame are analytically 
continued to the timelike region by changing q^^ to —q^ in the form factors. 

2. Zero-mode contribution 

In the nonvalence region P2 < < , the poles are at pi — Pion(^'^i) = [™i + — i^l/pf (from the 
struck quark propagator) and — p^^^-^{Ai) — [A^ + k^^ — ie]/pi (from the smeared quark-photon vertex), which 
are located in the upper half of the complex fc^-plane. In order to estimate the zero-mode contribution, we define 
a = P2 / Pi = 1 ^ 1^ / Pi =1^/3 and then the region P2 < < Pi corresponds toa<l — a;<lor equivalently 
Q < X < (3. That is, the zero-mode contribution((7+ — > 0) to the hadronic matrix element is obtained from the 
/3 0(i.e. X ^ 0) limit for the integration of the longitudinal momentum x. The fact that S'+j in Eq. (|20p is regular 
in the x — > limit implies no zero-mode contribution to f+{q'^). However, as one can see from Eq. (I2ip . both S^^ 
and include the terms proportional to l/a;(i.e. Pi), which are singular as x — > 0. Those singular terms in the 
perpendicular current may be the source of zero-mode contribution to the hadronic matrix element in the nonvalence 
region. 

When we do the Cauchy integration over k~ to obtain the LF time-ordered diagrams, we use Eq. pT|) to avoid the 
complexity of treating double p^-poles. As mentioned above, the zero-mode contribution comes from the p~{i = 1, 2) 
factors in S^^ = S;^^ + S^g^^. For instance, we define the zero-mode contribution to the l/lNqNj^-^NA^) term in Eq. (fTTj) 
having pj" = p~^^^(Ai) pole as 



[MUh.u.=^N]irn| 7^ • (28) 



d'^k 5,4 (Pi =Pio„( Ai)) 
ff^hj,^^ (27r)4 NqNA.NA, 

The zero-mode contributions to the other three terms in Eq. ([TT|l can be defined the same way as in Eq. (|28p to give 
the net zero-mode contribution TW^^. = I-^AiA^ (Pro„(Ai))]z.M. - [Mi^„^^{p^^^{Ai))]z.M. - [-A^iiA^ bron("^i))]z.M. + 
[■'^mim2 (?'ron("^i))]z.M. • Essentially, the nonvanishing zero-mode contributions in Eq. (j28p are summarized as follows: 



(1) 


Z = 


Pi > 


P2 , k in : 






lim 


/ 


d-^k Z 


i 




J nv 


{2n)^ NqNA.NA, 


167r3 


(2) 


Z = 


Plon ^nv '• 






lim 

/3-.0 


j 


d^k Z 


i 




J nv 


{2n)^ NqN A, N A, 


167r3 


(3) 


Z = 








lim 

/3-.0 


/ 


d^k Z 


i 




J nv 


{2n)^ NqNA.NA, 


167r3 



1 r a2 

dz / d^k 1 



1 

dz I d'k 



A2 


[zAl^ + (1 










A^ 


[zAj^ + (1 






m2j_z/ [z - 


-1) 


A^ 


[zAl^ + {l 





where the variable change x — fiz was made and r7i^_|_ = mf + p^^ and A^^ = A| + We should note that 

Z = p~,^k'^ PjonPii''- ~ 1' 2) in and other terms such as Z = PionPt ^^'^ P2onPt 8'° to zero in the /3 — > limit. 
Finally, we get the following nonvanishing zero-mode contribution to the trace term S^^^ in Eq. 



^^M. = lim 5i = 2pr(pii +P2±) = 2pr(2pii - Qj.), (30) 
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which in fact is common to other three terms in Eq. ((TT|) . After a httle manipulation, we finally get the following 
nonvanishing zero-mode contribution to the form factor f-(q^) in Eq. (jlSp as follows 

fZ.M.,2. _ , N f d^k ■^± f 1 I 1 

J- W J — Sa 2 9\/A 2 9\ / /o_^4 „2 AT i\r Ar at at at at at at 



{Ai^ - mi^)iA2^ - m2^) J {2n)^ qi V^g^Ai^A2 N^Na^N^ NgNiNA, NgNiN2 



87r2(A/-mi2)(A2^-TO2')yo \Ba,a,B. 



I dz{l - 2z) In 

^0 



(31) 



where 





= z{l- 




Hi 


-z)A?4 


-zAl 


BAim2 


^ z{l- 


-)qi- 


HI 


-z)A?-t 


' zmj, 


BrniA2 


= z{l- 


-)qi- 


HI 


— z)rn^ - 


^zAl 


Bmim2 


= z{l- 


^)qi- 


HI 


— z)mi - 


- zm2 



(32) 

Therefore, we get the LF covariant weak form factors in the (7+ = frame as /+^^°^(<z^) = /""'(g^) and jLFCov^^2-j ^ 

/™'(<Z2)+/Z.M.(^2)^ 



3. Effective inclusion of the zero-mode in the valence region 



In this exactly solvable covariant BS model, we find that while the matrix element of the plus current is exactly 
on-mass shell physical amplitude, that of the perpendicular current is the off-mass shell amplitude. As shown in our 
previous work [2q . we can relate the non- wave- function vertex to the ordinary valence wave function in the > 
frame using the iteration of the irreducible kernel involved in the bound state equation. In the frame, the 

nonvalence contribution in the g"*" > frame corresponds to the zero-mode contribution in the g"*" = frame. Thus, 
we can identify the zero-mode operator that is convoluted with the initial and final state valence wave functions 
to generate the zero-mode contribution. Our method can also be realized effectively by the method presente d by 
Jaus |3l] using the orientation of the light-front plane characterized by the invariant equation ui ■ x = [sgI l37j |. 
where oj is an arbitrary light-like four vector. The special case uj = (1, 0, 0, —1) corresponds to the light-front or null 
plane u; ■ x = =0. While the exact on-shell amplitudes (such as A^+) should not depend on the orientation of 
the light-front plane, the off-shell matrix elements(such as Ai'^) acquire a spurious u; dependence. This problem is 
closely associated with the violation of rotational invariance in the computation of the matrix element of a one-body 
current. In order to treat the complete Lorentz structure of a hadronic matrix element the authors in [sil [36j have 
developed a method to identify and separate spurious contributions and to determine the physical, i.e. co independent 
contributions to the hadronic form factors. Below, we summarize the result of zero-mode contribution obtained from 
the method by Jaus [Slj] and discuss the equivalence with our result of zero-mode contribution. 

By adopting the oj dependent light-front covariant approach as in [3l|, [s^ , we derive the light-front covariant form 
of the form factor /-(g^), which effectively includes the zero- mode contribution in the valence region. In order to do 
this, we first decompose the four vector in terms of P = (Pi + P2), g, and uj with uj = (1, 0, 0,-1) as follows [31|] 



= pf'A[^^ + (7^4'^ + ^T^^^'Cf ^ • (33) 



UJ ■ P 



The coefficients in Eq. ((33)) are given by 





UJ 


Pi 




UJ 


P 




1 


{pi 




- ? 






= Pi 


■P 



X 

2' 



[q-P)- 



uj- P J 2 g2 ' 

Pi-P- P'AY' - q ■ PA'^2^ =Z2- Nq, (34) 



where Nq = fc^ — m? and 



Z2 = x{Ml - Ml) +ml- ml + {1 - 2x)M^ - [q^ + q ■ P]^^^. (35) 
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Note that only the coefhcient C^"'^^ which is combined with lu^ depends on p-^ (i.e. zero-mode). In this exactly solvable 
BS model, the zero-mode contribution from is exactly opposite to that from Nq, i.e. 



Z.M. 



Pi 



N 

= -I[Nq]z.M.. 



dz In 



(36) 



Furthermore, the zero-mode contribution I[Nq]z.M. from Nq is exactly the same as the valence contribution I[Z2 



val 



from Z2, where I[Z2]vai is given by 

1 dx f 

I[Nq]z.M. = I[Z2]vai = -77-3 / -j / d'^^±Xi{x,k±)x2{x,k'j^)Z2. 



(37) 



From the identities in Eqs. (|36p and ([37]), the replacement Nq — > Z2{ov equivalently —Z2) in the spurious 

Lu dependent (i.e. the zero-mode related) term C^^'' in Eq. (j34p becomes covariant, i.e. free from any to dependence. 
Effectively, the zero-mode contribution from pi in the valence region can be given by Eq. (j37|) . Using this, we can 
effectively include the zero-mode contribution from the second term Pj^qj^ in Eq. (|30|) in the valence region. On the 
other hand, since the first term PiPi± in Eq. (|30p has a tensor structure, we need the tensor decomposition (3l| 



1 



P^T^A^^^ + {P^'q" + q^P'')A'^^^ + q^'q'^Af'^ 



1 



UJ ■ P 



(P^'cj" +uj^'P'')B[ 



(2) 



LU ■ P 



iq^iu'' +LU^q'')C[ 



(2) 



(lo-P) 



(38) 



where 



A 



(2) 



(k± • q±) 



A. 



(2) 



(l)l2 



A 



(2) 



.(l).(l) .(2) 
^1 ^2 ' ^4 



SI 



(2) 



^(1)^(1) 4(2) ^(2) 



4(1)^(1) , 9- 4(2) ^(2) 



(1)12 



p2 



^2 A ' 

(g ■ P? 

q2 



A\ 



(2) 



(39) 



We note that the coefficients ' {i = 1, • • • ,4) are related with = + or _L components and B[ ' with z/) = 
(+, — ). According to our power counting rules mentioned above, those terms are zero-mode free. On the other hand, 

(2) (2) 

the coefficients and are related with {iJ.,i') = (— ,-L) and (—,—), respectively. That is, the C terms are 

(2) (2) 

related with the zero- mode contributions. Specifically, and are related with the zero-mode contribution to 
the perpendicular and minus components of the currents, respectively. The zero-mode contribution from PiPii_ is 

(2) 

thus related with ' and the effective inclusion of the zero-mode in the valence region can be achieved by setting 



0. This leads to 



A«Z2 



-^A\ or pij 



Ai}^Z2 



1 ■ Pa^) 

„2 ^1 



(40) 



In summary, the zero-mode contribution from Sz m given by Eq. (j30p can be expressed in terms of the zero- mode 
operator convoluted with the initial and final state LF vertex functions: 



US' 



1 



\Z.M. 



dx 



167r3 Jo (1 - 



d^k±xi{x,k±)x2{x,k'j_)\ -4q 



Ai'^Z2 



q-P 



(2) 



2q±Z2 



(41) 



as expected from our effective method presented in our previous work 26]. Consequently, the LF covariant form of 
the form factor f-{q^) is obtained as 



j (^ilrj.-^ J d'^'^±Xi{x,ki_)x2{x,k' ±)^^-x{l -x)MI - k\ - miniq + (m2 - ■mq)Ai 



, (k± • qj 



k± ■ q± 



[M| ^ {l-x){q^ + q- P) + 2xM^ 



-(1 - 2x)Ml - 2(mi -mg)(mi +7712)] >, 



(42) 
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' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' 

B --> n (in fermion field theory model) 




5 10 15 20 25 30 

FIG. 1: The normalized weak form factors f±{q'^)/ f+{0) for 
B ^ n semileptonic decays obtained from the exactly solvable 
covariant BS model of fermion field theory. 



Ml 



where q ■ P 

In Fig. [Tl we show our results of the normalized weak 
form factors f±{q^)/f+{0) for the semileptonic B ^ tt 
decay obtained from the exactly solvable covariant BS 
model of fermion field theory. The used model parame- 
ters for B and tt mesons are Mb — 5.28 GeV, AIt^ — 0.14 
GeV, rrib = 4.9 GeV, TO„(d) = 0.43 GeV, = 10 GeV, 
and A„(j;) = 1.5 GeV. These parameters are fixed from 
the normalization conditions of the tt and B clastic form 
factors at g2 = [Hi. 

The solid line represents the form factors /+(g^)//+(0) 
obtained from the manifestly covariant result [Eq. (fT3H 



and from the full LF calculation /^'"^"^(g^) = /^'(g^) 
[Eq. (P7| ] in the g+ = frame. Since the two results are 
in complete agreement with each other, we depict them 
by the single solid line. The dashed line represents the 
form factors — /_((7^)//+(0) obtained from the manifestly 
covariant result [Eq. (|13p ] and from the full LF calcula 



pends on the orientation of the light-front plane or equiv- 
alently receives zero-mode contributions. (2) Our power 
counting rule for is very efficient method in identifying 
the zero-mode terms such as Pi and pj~Pi_L that appear 
in the perpendicular currents and PiPi appearing in the 
minus current. (3) We explicitly show that the uj depen- 
dent(i.e. zero-mode related) coefficients c[^\c[^\ and 

Cj^-* correspond to pj~,pj~Pi_L and PiPi terms, respec- 
tively. These features in our approach should be distin- 
guished from the approach presented in [^ij . 

While the manifestly covariant BS model of fermion 
field theory model is good for the qualitative analysis of 
semileptonic decays, it is still semi-realistic. We thus dis- 
cuss more phenomenological LFQM and the LF covariant 
form factors within our LFQM in the following sections. 



III. MODEL DESCRIPTION 

The key idea in our LFQM [2^ for mesons is to 
treat the radial wave function as a trial function for 
the variational principle to the QCD-motivated effective 
Hamiltonian saturating the Fock state expansion by the 
constituent quark and antiquark. The QCD-motivated 
Hamiltonian for a description of the ground state meson 
mass spectra is given by 



fc2 



[H0 + Vgg]\'^ 



fc2 + Vqq 



I nlml ' 



(43) 



, fcz) is the three- momentum of the con- 



where k = (kj 

stituent quark, Mqq is the mass of the meson, and Y^'^^ 
is the meson wave function. In this work, we use two 
interaction potentials Vqq] (1) Coulomb plus harmonic 
oscillator(HO) and (2) Coulomb plus linear confining po- 
tentials. The hyperfine interaction essential to distin- 
guish pseudoscalar(0 ^) and vector(l^^) mesons is also 
included; viz., 



Vqq Fo + Vl,. 



Vc, 



4a^ 
3r 



2S„ -S 



3 mqrriq 



tion /. 



LFCov 



Z.M, 



(44) 

where Vconf = hr{r^) for the linear (HO) potential and 
lq^)\EQ. raiTnThe ' ^^i] " V4(-3/4) for the vector (pseudoscalar) me- 



— frame. Here again, the two results are in complete 
agreement with each other. The dotted line represents 
only the valence contribution to -~f-{q^)/f+{fi). The 
difference between the dashed and dotted lines amounts 
to the zero- mode contribution to the form factor f-{q^). 

Although our result for the /LFC ov^ 21 jg essentially the 
same as that obtained from Jaus [31], the distinguished 
features of our approach in deriving the LF covariant 
form factor may be summarized as follows: (1) We sepa- 
rate the trace term into the on-shell propagating part 
and the instantaneous part S^^^^^^ which enable us to clas- 
sify the on-shell and off-shell matrix elements explicitly. 
From this one can easily find which matrix element de- 



son. Using this Hamiltonian, we analyze the meson mass 
spectra and various wave-function-related observables, 
such as decay constants, electromagnetic form factors of 
mesons in a spacelike region, and the weak form factors 
for the exclusive semileptonic and rare decays ofpseu- 
doscalar mesons in the timelike region [23, [H, [2^, [l^, [H, 

The momentum-space light-front wave function of the 
ground state pseudoscalar and vector mesons is given by 



100 



f (a;i,ki_L, Ai) ^Ti'l'll^{xi,'ki±)(l){xi,'ki±), (45) 



where (t){xi,'kii_) is the radial wave function and ^AiAs 
is the spin-orbit wave function that is obtained by the 
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interaction- independent Melosh transformation from the 
ordinary spin-orbit wave function assigned by the quan- 
tum numbers J^'-^ . The model wave function in Eq. (|45|) 
is represented by the Lorentz-invariant internal vari- 
ables, Xi — pI / , = Pi± — ±_ and Ai, where 
= (p+,p-,p^) = (po + p3,(A/2 + p2_)/p+,p_L) is 
the momentum of the meson A/, and and are the 
momenta and the helicities of constituent quarks, respec- 
tively. 

The covariant forms of the spin-orbit wave functions 
for pseudoscalar and vector mesons are given by 



above procedure, we find an analytic form of the mass 
eigenvalue given by 




as 



2/32 

8/3 , 32/33(S, -S^ 



(50) 



7^' 



-'"Ai(Pl)75'yA2(P2) 



Ai A2 



V2M0 



n 



-"Aibi) 



e-(pi-p2) 



V2M0 



(46) 



where Mq 



m~)/xi is the boost invariant meson mass square ob- 
tained from the free energies of the constituents in 
mesons, and e^{Jz) is the polarization vector of the vector 
meson [38|. The spin-orbit wave functions satisfy the re- 
lation X^AiAa '^\i\l'^\iX2 ~ ^ both pseudoscalar and 
vector mesons. For the radial wave function 0, we use 
the same Gaussian wave function for both pseudoscalar 
and vector mesons: 



02 



^exp(-fcV2/32), (47) 



where /3 is the variational parameter. When the longi- 
tudinal component is defined by — {x — 1/2) Afg + 
(m2 — m2)/2Afo, the Jacobian of the variable transfor- 
mation {x, kx} —tk= (k^, kz) is given by 



dx 



Mq 
4xiX2 





r 2 21 






m'l — r?i2 




[ M§ \ 





(48) 



Note that the free kinetic part of the Hamiltonian Hq = 
is equal to the free mass operator 
Mq in the light-front formalism. 

The normalization factor in Eq. (I47|) is obtained from 
the following normalization of the total wave function: 



dx 



1. 



(49) 



where 1 = 



and Ki{x) is the modified Bessel func- 



We apply our variational principle to the QCD-motivated 
effective Hamiltonian first to evaluate the expectation 
value of the central Hamiltonian Hq + Vo, i.e., {(j)\{Ho + 
Vo)\4>) with a trial function (j){xi,\ii±) that depends 
on the variational parameter /3. Once the model pa- 
rameters are fixed by minimizing the expectation value 
{(f)\{HQ + Vb)|0), then the mass eigenvalue of each meson 
is obtained as Mqq = {(j)\{Ho + Vqq)\(f)). Following the The decay constants of pseudoscalar and vector mesons 



tion of the second kind. The upper and lower components 
of the column vector in Eq. (|50p represent the results for 
the linear and HO potential models, respectively. Mini- 
mizing energies with respect to /3 and searching for a fit to 
the observed ground state meson spectra, our central po- 
tential Vq obtained from our optimized potential parame- 
ters (a = -0.72 GeV, b = 0.18 GeV^, and as = 0.31) Q 
for Coulomb plus linear potential was found to be quite 
comparable with the quark potential model suggested by 
Scora and Isgur [39| where they obtained a = —0.81 GeV, 
b = 0.18 GeV^, and = 0.3 - 0.6 for the Coulomb plus 
linear confining potential. A more detailed procedure for 
determining the model parameters of light- and heavy- 
quark sectors can be found in our previous works ^24 . i25| . 
In this work, we obtain the new variational parameter 
f3cb for the bottom-charm sector and predict the mass 
eigenvalues of the low-lying Be and B* states. Our new 
prediction of Mb^ — 6459 [6351] MeV obtained from the 
linear [HO] potential model is in agreement with the data, 
A^gP = (6276 ± 4) MeV ^ within 3% error. We also 
predict the unmeasured mass of B* as Mb* = 6494 [6496] 
MeV for the linear [HO] potential model. Although it is 
generally believed that the linear potential is preferred 
between quark and antiquark in the heavy meson sys- 
tem, our result of the spectrum computation indicates 
that the HO potential is also viable and thus leads to the 
further investigation. We use both HO and linear poten- 
tials to compute the decay processes presented below and 
identify the physical observables sensitive to the choice 
of potential. 

Our model parameters {mq,f3qq) and the predictions 
of the ground state meson mass spectra obtained from 
the linear and HO potential models are summarized in 
Table U and in FigJ5J respectively, compared with the 
experimental data [40| • Our prediction of the r]b meson 
obtained from the linear [HO] potential model, M^^ = 
9657 [9295] MeV slightly overestimates [underestimates] 
the very recent data from the Babar experiment, Af^^P = 
9388.9t^;3(stat) ± 2.7(syst) MeV [4l|. OveraU, however, 
our LFQM predictions of the ground state meson mass 
spectra are in agreement with the data px| within 6% 
error. 
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TABLE I: The constituent quark mass [GeV] and the Gaussian parameters /3 [GeV] for the hnear and HO potentials obtained 
by the variational principle, q = u and d. 



Model 


niq 


nis 


rric 


TUb 






(5ss 






Pec 




Psb 




Pbb 


Linear 


0.22 


0.45 


1.8 


5.2 


0.3659 


0.3886 


0.4128 


0.4679 


0.5016 


0.6509 


0.5266 


0.5712 


0.8068 


1.1452 


HO 


0.25 


0.48 


1.8 


5.2 


0.3194 


0.3419 


0.3681 


0.4216 


0.4686 


0.6998 


0.4960 


0.5740 


1.0350 


1.8025 



(9295) 


Tl^(9389) 


(9657) 


(9558) 


Y(9460) 


(9691) 


(6351) 


B.(6276) 


(6459) 


(6496) 


B»^(?) 


(6494) 












(5378) 
(5235) 


B^(536S) 
B(5279) 


(5375) 
^(5235) 


(5471) 
(5349) 


B«^(5413) 
B«(5325) 


(5424) 
(5315) 


(3128) 


tl_,(2980) 


^:3171) 


(3257) 


J/i/(3097) 


(3225) 


(2005) 
(1821) 


DJ1968) 
D(1869) 


(2011) 
(Kid) 


(2150) 
(2024) 


D*,(2112) 
D«(2010) 


(2109) 
(1998) 


(470) = 


Tl'(958) 
11(547) 
K(494) 
^(MO) 


= (478) 


(875)= 

k'(892) 


(K1020) 

P(~ 

m<782) 




(HO) 


(Exp.) 


(Linear) 


(HO) 


(Exp.) 


(Linear) 



FIG. 2: (color online). Fit of the ground state meson 
masses[MeV] with the parameters given in TableU] The (p, tt), 
(77,77'), and {Lu,(f>) masses are our input data. The masses of 
[ijj — (j)) and (77 — 77') were used to determine the mixing angles 
oi u) — (f> and rj — r]' [3], respectively. 



are defined by 



\qrq\V{P,h)) 



fvMve^{h). 



(51) 



In the above definitions for the decay constants, the ex- 
perimental values of pion and rho meson decay constants 
are « 131 MeV from vr ^ ^i/ and fp w 220 MeV from 
p — > e+e^. 

Using the plus component (/i = +) of the currents, 
one can easily calculate the decay constants. The explicit 
forms of pseudoscalar and vector meson decay constants 



TABLE II: Bottom-charm meson decay constants(in unit of 
MeV) obtained from the linear [HO] parameters. 



Linear [HO] 


f5] f8] 


[42] 


[43] 


[44] 


[45] 


fB, 377[508] 


360 433 


500 


460 ± 60 


517 


410 ± 40 


/s- 398[551] 


- 503 


500 


460 ± 60 


517 





are given by 
fp 
fv 



2V6 
2V6 



dx (i^kj_ A 
167r3 

dx d^k 1 



(a;,kx), 



167r3 



A 



2ki 



,(52) 



where A = X2'mi -f Xim2 and A^o = -^-^^o + "^tii + '7i2- 
Here, only Lz = Sz = component of the wave func- 
tion contributes. We note that the vector meson decay 
constant fv is extracted from the longitudinal {h — 0) 
polarization. 

In Table |lTl we present our predictions for the decay 
constants of fs and fs- and compare with other model 
calculations [1,11, [H, 0, Hi, SI ■ The decay constants 
for other light- and heavy-mesons have been predicted 
in our previous works [2J, [28l . [38j and found to be in 
good agreement with experimental data. While our pre- 
dictions of the decay (constants for the light-light and 
heavy-light systems p3 . [28l Isst are not sensitive to the 
choice of potential(linear or HO), the decay constants for 
heavy-heavy systems such as {Be, B*) and {r/h, T) in [llj 
are quite sensitive to the choice of potential. Thus, the 
experimental measurements for the decay constants of 
(Be, B*) and (77^, T) mesons may distinguish between the 
linear and HO potentials within our LFQM. 

The process-independent quark distribution ampli- 
tude (DA) (f)p(y){x) for pseudoscalar (vector) meson is 
the probability amplitude for finding the qq pair in the 



meson with Xa = x and Xs 



1 



to our LF valence wave function 



It is directly related 



IGtt- 



-^'(a;,kj 



(53) 



The k^ integration in Eq. (|53|) is cut off by the ultraviolet 
cutoff A implicit in the wave function. The dependence 
on the scale A is then given by the QCD evolution equa- 
tion [ill and can be calculated perturbatively. However, 
the DAs at a certain low scale can be obtained by the nec- 
essary nonperturbative input from LFQM. Moreover, the 
presence of the damping Gaussian factor in our LFQM 
allows us to perform the integral up to infinity without 
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loss of accuracy. The quark DAs for pseudoscalar and 
vector mesons are constrained by 

''\p^y^{x)dx = ^. (54) 

We show m Fig. [3] the normaUzed quark DAs ^{x) = 
{2y/6/fp)(j){x) for D (dotted Une), B (dashed hne), 
(dot-dashed hne), and Be (sohd hne) mesons obtained 
from the hnear (upper panel) and HO (lower panel) po- 
tential parameters, respectively. In Fig. |4l we also show 
the normalized quark DAs for rjc (thin lines) and ryt (thick 
lines) mesons obtained from the linear(solid lines) and 
HO (dashed lines) potential parameters. While the two 
model predictions for the heavy-light systems such as 
(D, B, and Bg) are not much different from each other, 
the HO potential model predictions for the heavy-heavy 
systems such as {rjc, Be and rjb) give somewhat broader 
shapes than the linear potential model predictions. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



IV. SEMILEPTONIC DECAYS OF THE 
MESON 

The relevant quark momentum variables for P{qiq) 
P{q2q') transitions in the q^ = Q frames are given by 

p+ = XiP^ , = X2P^, Pl± = XiPi± - k^, 

P2± = XiP2± - k'L, 



P2 = XlPl, Pg,^X2Pi 



Pq± = X2P1 



k' 



(55) 



k±, Pq'± = X2P2A 

where 2:1=2: and 2:2 = 1 — a; and spectator quark requires 
that pq — and p^^ = Pq'±- Taking a Lorentz frame 
where Pij^ = and P2± = amounts to k'^ = 

kj_ (1 - x)q±. 

The hadronic matrix element of the plus current 
M+ = {P2\V+\Pi) in Eq. (g]) is then obtained by the 
convolution formula of the initial and final state LF wave 
functions in the valence region: 



M+ = 



1 '• d^k^ 



dx 



167r3 



(/)2(a;, \i'±_)4>i{x, kj_) 



J2 7^ooj^iA^^+^iA^^oo__(5g) 



Substituting the covariant form of the spin-orbit wave 
function for pseudoscalar meson given by Eq. ()46|) into 
Eq. (iniD yields 



7W+ = - 



dx 



(x.k' 



1(2;, kj_) 



/o J 1671-^ 

X Tr[7^(^2 + ^2)7^ 



2xPj^MoM'a 



where 



M'n = 



(mi 



m?)2 = 







2:(1 


-x)' 


/k'i 





x{\ — x) 



(58) 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 08 0.9 1 



FIG. 3: The normalized distribution amplitudes for D, B, 
Bs, and Be mesons obtained from the linear (upper panel) 
and HO (lower panel) potential parameters. 



and Ai = (1 — x)mi + xrriq (i = 1, 2). After some manip- 
ulation, the trace term in Eq. (I57p is reduced to 



Tr[7^(|*2 + ^2)7^(1*1 + ™i)7^(^9 - ™)] 

^(ki.kl+^1^2). 



(1-.T) 



(59) 



Finally, the form factor f+{q^) obtained from the valence 
contribution in the q'^ — Q frame is given by 

r ( 2-. _ f\ f d^^i- Mx,'i^±) 02 (a:, k'^) 

X {AiA2 + 'k±-'k'±). (60) 
We should note in the trace calculation of Eq. ([55)1 that 
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FIG. 4: The normalized distribution amplitudes for -qc (thin 
lines) and rjb (thick lines) mesons obtained from the linear 
(solid lines) and HO (dashed lines) potential parameters. 



the internal momenta of the valence quarks carried inside 
mesons are all on- mass-shell {pf = ml). Nevertheless, 
the LF valence contribution to the form factor f+iq^) is 
shown to be equivalent to the covariant result as shown 
in Sec. II. 

Comparing the manifestly covariant form factor f+{q^) 
in Eq. ^ and our LFQM result f+{q^) in Eq. jM]), we 
find the following relations for the LF vertex functions 
between the two models: 



'2N 



/2N 



Xi{x,'kA_) _ 01 (a;, k^) 
1 - X 

X2(a;,k'i) 



1 



4>2{x,'k' j_) 



(61) 



We should note that the zero-mode operator included in 
Eq. ([42]) is independent from the choice of radial wave 
function. 

Applying the relation in Eq. dH]) to Eqs. ^ and (g^, 
we get the following LF valence contribution /^'(q^) and 

the LF covariant solution fly'^\q^) including both the va- 
lence and the zero-mode contributions within our LFQM: 



r-\q^) = I {l-x)dx 



d^ki (/)i(a;,k^) (j)2{x,'k' i_) 



J2 



'^^ \Ml + Ml - 2(mi - m,-)(m2 - m,-)] j>, 



-(1 — x)M'1 + (to2 — mq)Ai — mq{mi — mq) 



(1 - 2x)Mf - 2(mi - mg)(mi + ma)] }, 



x{l — x)Mi — k^ — mirrig + (rria — mq)Ai 



k± qj 



[Ml -{l-x){q^ +q- P) + 2xMl 



(62) 



where Mi and M2 are the physical masses of the initial 
and final mesons, respectively. 

V. RADIATIVE B* ^ B^-y DECAY 

In addition to semileptonic decays, the radiative decays 
of vector mesons can be anlayzed within our LFQM [2J, 
[28j . In this work, we thus calculate the decay rate for 
B* Bc'^ transition. 

In our LFQM calculation of B* — > i3c7 process, we 
first analyze the virtual photon (7*) decay process, cal- 
culating the momentum dependent transition form fac- 
tor, Fs^Baiq^)- The lowest-order Feynman diagram for 
V P7* process is shown in Fig. [5] where the decay from 



vector meson to pseudoscalar meson and virtual photon 
state is mediated by a quark loop with flavors of con- 
stituent mass mi and niq. 

The transition form factor Fb-bS^^) f^^' ^(-^i) ~^ 
P{P2) + l*{q) is defined as [H 

{P[P2W\V{Pi,h)) = iee'^''P-e,{Pi,h)qpPi^Fvp{q^), 

(63) 

where the antisymmetric tensor e^"''"' assures electro- 
magnetic gauge invariance, g = Pi — P2 is the four- 
momentum of the virtual photon, and e^{Pi,h) is the 
polarization vector of the initial meson with the four- 
momentum Pi and the helicity h. The kinematically al- 
lowed q^ (momentum transfer squared) ranges from to 
9max = {Mv - MpY. The decay form factor Fb'bM'^) 
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FIG. 5: The lowest order Feynman diagram for V 
process. 
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can also be obtained in the ~ Q frame with the trans- 
verse {h = ±1) polarization and the "+" -component of 
the currents without encountering zero-mode contribu- 
tions [S^l and then analytically continued from the space- 
like region where the form factor is given by ^^.^^(q^^) 
to the timelike > region by changing q'^ to — in 
the form factor. 

The hadronic matrix element of the plus current 
{P{P2)\V+\V{Px,h)) in Eq. ^ is then ob- 
tained by the convolution formula of the initial and final 
state LF wave functions in the valence region: 



M+p 



M 



VP 



3 



dx 



167r3' 



2{x,k' ±)(l)i{x,k±) 



"A, A 



^K.' (64) 



where ecj is the electrical charge for jth quark flavor. 
Substituting the covariant forms of the spin-orbit wave 
functions for pseudoscalar and vector mesons given by 
Eq. ([^5]) into Eq. and comparing it with the right- 
hand side of Eq. i.e. eP^ Fvp{q^)q^' / \/2 where 



R _ 



(Ix + i<lyi we could extract the one-loop integral, 
(Pk^ 0i(a;,k_L) 02(a;,k'_L) 



/(mi,TOg, given by 
I{mi,m,q,q^)= J dx J 



8^3 + 

(k 



il-x)iM + —[k: 



\2 



q- 



(65) 



where A4o = Mq -|- mi -|- nig. The decay form factor 
Fb-bAi^) is then obtained as [11] 



B'B, 



(q^) = ei/(mi, m^, q"^) + e2/(m^, mi,q^). (66) 



The coupling constant gs^B^-y for I'cal photon (7) case 
can then be determined in the limit q^ — )■ 0, i.e., gs^B^-y — 
0) . The decay width for V Pj is given by 



Fb'bAq^ 



yBiBai'^'y' 



(67) 
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FIG. 6: (color online). The weak form factors f+{q^) (solid 
line) and fo{(f) (dashed line) for Be ^ D (upper panel) and 
Be gc (lower panel) semileptonic decays obtained from 
the linear (black line) and HO (blue line) potential parame- 
ters. The circles represent the valence contributions /o'''(g^) 
to Mq^). 



where a is the fine-structure constant and 



Af^ )/2Mb* is the kinematically allowed energy of the 
outgoing photon. 



VI. NUMERICAL RESULTS 

In our numerical calculations of exclusive Be decays, 
we use two sets of model parameters (m, (3) for the linear 
and HO confining potentials given in Table [T] obtained 
from the calculation of the mass spectra. Although our 
predictions of ground state heavy meson masses are over- 
all in good agreement with the experimental values, we 
use the experimental meson masses [13] in the computa- 
tions of the decay widths to reduce possible theoretical 
uncertainties. We also use the central values of the CKM 
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FIG. 7: (color online). The weak form factors f+{q'^) and 
/o(g^) for Be B (upper panel) and Be Bs (lower panel) 
semileptonic decays obtained from the linear and HO poten- 
tial parameters. The same line codes are used as in Fig. [6] 




ry-'/Gt-V-'y 

FIG. 8: (color online). The weak form factors f+{q^) and 
/o(q^) for rjt Be semileptonic decay obtained from the lin- 
ear and HO potential parameters. The same line codes are 
used as in Fig. |S1 



matrix elements, 

\Vub\ = 0.00393, IKhl =0.0412, 

\Vcd\ = 0.230, 114.1 = 1.04, (68) 

quoted by the Particle Data Group (PDG) 

In Figs. [H] and [7] we show the g^-dependence of the 
LF covariant weak form factors /+(g^) (solid lines) and 
/o((7^)(dashed lines) in the whole kinematical ranges 
for the CKM-suppressed (enhanced) semileptonic Be 
D{ric) (Fig. El) and B{B,) {Vig.^ decays obtained 

from both linear(black lines) and HO (blue lines) poten- 
tial models. The circles represent the valence contribu- 
tions /^'(g^) to fo{q'^)- That is, the difference between 
fo{q^) and /™'(<7^) represents the zero-mode contribu- 
tion to /o(g^). 

The kinematical ranges for Be —> D{r]c) decays induced 
by 6 ^ u{c) transitions with the c quark being a specta- 
tor are considerably broader than those for Be B{Bs) 
decays induced by c — *■ d{s) transitions with the b quark 
being a spectator. The form factors f+iq^) and fo{q^) at 
the zero-recoil point (i.e., q^ — (Zmax) correspond to the 
overlap integral of the initial and final state meson wave 
functions. The maximum-recoil point (i.e., q^ = 0) corre- 
sponds to a final state meson recoiling with the maximum 
three-momentum \Pf\ = {Mg — M'j) /2Mb^ in the rest 
frame of the Be meson. Especially for the Be ^ D decay, 
the light u quark in D meson will typically recoil with 
the momentum comparable to or larger than the c quark 
mass due to the large recoil effect for Be ^ D decay. In 
order for the final D meson to be bound, there must be 
a correspondingly large momentum transfer to the spec- 
tator c quark. Thus, the overlap between the initial and 
final meson wave functions at the maximum-recoil point 
is limited and yields smaller value of /-i-(0) for Be ^ D 
decay than that for other processes. We also note that 
one cannot apply the heavy quark symmetry to the sys- 
tem with the two heavy quarks, due to the flavor sym- 
metry breaking by the kinetic energy terms as discussed 
in [47|. As for the zero-mode contributions, we find that 
the zero- mode contributions to /o((7^)(or f-{q^)) for the 
Be — > D,Bc B and Be Bg processes are rela- 
tively larger than that for the Be rjc process. For 
the Be — > rje transition, the zero-mode contributions to 
fo{q^) obtained from both linear and HO potential mod- 
els are almost suppressed in the whole kinematical range 
and moreover the values of /o(9max) almost converge to 
a single value regardless of the choice of potential. 

Although there already exist various model predictions 
on the above Be semileptonic decays, the predictions of 
the semileptonic 77^ — s- Be decay is not reported yet as far 
as we know. We thus show in Fig.[8]the g^-dependence of 
the weak form factors f+{q^) and foiq'^) for the semilep- 
tonic 77b — > Be decay obtained from both linear and HO 
potential models. The same line codes presented in Fig. [HI 
are used in Fig. [51 While the linear and HO poten- 
tial models give similar decay constants for heavy-light 
mesons {D,B,Bs) and rje meson [11], they predict quite 
different values of Be and rji,, e.g., /^^ = 507 MeV and 
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FIG. 9: DifTerential decay widths {\ l\VQ^Q^f)dV j dq^)(\-a 
units of 10"^^ GeV"^) for Be ^ D£ue and Be r]dvi pro- 
cesses obtained from the linear and HO potential parameters. 



FIG. 10: Differential decay widths {l/\VQ^Q^\^)dV /dq^){m 
units of 10"^^ GeV"^) for Blvi and Be ~* Bsln pro- 

cesses obtained from the linear and HO potential parameters. 



897 MeV for the linear and HO potential models [2^ , re- 
spectively. This results in sizable differences between the 
two models for the predictions of and /o('Z^) in 

the rjb ^ Be decay. Since the linear potential model pre- 
diction of the quark DA for rjb is narrower than the HO 
model prediction (see Fig. 5]), the overlap between the 
initial and final meson wave functions at the maximum- 
recoil point (i.e., = 0) produces smaller values of 
f+{— fo) for the linear potential model than for the HO 
potential model. The experimental measurement of this 
process may also distinguish between the linear and HO 
potential models within our LFQM. The zero-mode con- 
tribution to /_(g^)(or foiq'^ j) is again quite suppressed 
in the whole kinematical range as in the case of Be — > rjc 
process. 

In Figs. [SlfTTl we show the differential decay widths 
dr/dq^ for the D{r]c)ivi (Fig. [9]), Be B{Bs)lvi 

(Fig. [ini) and ry?, — > Bcivt (Fig. [Tl|) processes obtained 
from the linear and HO potential parameters. The line 
codes are described in each figure. We should note that 



/ = T (I^inear) - 
/ = c (H(i) 
/ = T (HO) 



FIG. 11: Differential decay width (I/IVq^q^ |2)dr/dg2)(in 
units of 10~^^ GeV~^) for r\h Bc£ue process obtained from 
the linear and HO potential parameters. 
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the minimum value of the form factor depends on the 
actual final lepton and it is given (neglecting neutrino 
masses) by the lepton mass as q^^^ — mj. Although the 
difference between the linear and HO model predictions 
are not very large for the Be {B,Bs) processes, they 
are quite different for other processes, especially for the 
r]h —> Be process. Since the constituent masses of b- and 
c quarks are common to both linear and HO potential 
models, the difference of the decay rates for the r]b Be 
process seems to come from the different choice of the 
variational f3 parameters. We note, however, that the 
difference of the decay rates between the two models are 
significantly reduced for the heavy r lepton case. 

In Table IIII[ we summarize our results for the weak 
form factors /+ and fo at = and q^^,^ and the de- 
cay widths Ti of the semileptonic Be {D, rje, B, Bs)ii'e 
and rji, Betvi{(. — e,/i,T) decays in comparison with 
other theoretical model predictions [1, [^, 

M m, M m, 

HgI . [itI . lisj . The subscript for the decay width Ti repre- 
sents the result for P Ptvi decay where the final lep- 
ton is ^ = e, /.t or r. For the decays induced by 6 — > u(c) 
transitions such as Be D, rjc and r]{, Be decays, 
we take Tg ~ T^ with the massless lepton limit since 
the muon mass effect is negligible for these transitions 
with large kinematic ranges. For the decays induced by 
c d{s) transitions such as Be B{Bs) decays, is 
about 5% smaller than Fe in our model predictions. For 
the Be ~* D decay, our predictions of the form factor 
at the maximum-recoil point are rather smaller than 
other quark model predictions. The upcoming experi- 
mental study planned at the Tevatron and at the LHC 
may distinguish these different model predictions. For 
the Be r/cB and Bg semileptonic decays, our pre- 
dictions are quite comparable with those of the quasipo- 
tential approach to the relativistic quark model d, [l3| , 
the relativistic quark- meson model [l5| , and the nonrela- 
tivistic quark model ^1^. It may be noted, however, that 
the predictions of the quark model based on an effective 
Lagrangian describing the coupling of hadrons to their 
constituent quarks fs'] as well as the covar iant LFQM [l^ 
are quite different from other model predictions including 
ours. 

Finally, in order to analyze the total rate for the ra- 
diative _B* — > i?c -t- 7 decay, the masses of the Be and 
B* mesons must be specified. Although we predicted the 
above two meson masses in Fig. [21 we use the the central 
value of the experimental data M'^'^ = 6.276 GeV [i^l 
to reduce the possible theoretical uncertainties. For the 
unmeasured B* meson mass, we take some range of the 
B* meson mass, i.e., 10 MeV < Am(= M^. - MbJ < 
220 MeV. The upper value of Am(i.e., M^. = 6496 
MeV) is chosen to be corresponding to our predictions, 
Mb* = 6494 MeV and 6492 MeV, obtained from the 
linear and HO potential models, respectively. 

In Fig. I12i we show the momentum-dependent form 
factor FB'sA'f) (upper panel) for the radiative B* 
Be^* decay and the dependence of r(_B* Be^) on Am 
(lower panel) obtained from the linear (solid line) and 



0.26 





lOO 150 
Am(=lVI - M )LIVleVJ 



FIG. 12: The transition form factor Fb*Bc{<1^) {upper panel) 
for — » -Bc7* and the dependence of r(_B* ^ _Bc7)(lower 
panel) on Am = Mb* — Mb^ obtained from the linear and 
HO potential parameters, where the central value of the ex- 
perimental data M'^'^ = 6.276 GeV 40] is used. 



HO potential (dashed line) parameters. For the transi- 
tion form factor Fb^bA^^)^ have performed the ana- 
lytic continuation of Fb^bAi^) from the spacelike region 



{q < 0) to the physical timelike < g 



< 



where 



^max 



{Mb- — Mb^Y represents the zero recoil point of 
the Be meson. The coupling constant qb-b^ is obtained 
at the g^ = point that corresponds to the Be meson 
recoiling with the maximum three-momentum in the rest 
frame of the B* meson. In our model calculation, the cou- 
pling constant itself is independent of the physical masses 
of the mesons and our prediction is gB*B^ = 0.273 [0.257] 
GeV^^ for the linear [HO] potential model. Our predic- 
tions are quite comparable with the result from the QCD 
sum rule approach [49|, g^f^^ = 0.270 ± 0.095 GeV-^. 
As one can see from the lower panel of Fig. [T^l the de- 
pendence of r(_B* Bej) on Am is quite sensitive to 
the mass of the B* meson, e.g., our linear [HO] poten- 
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TABLE III: Form factors /+ and fo evaluated at — and q, 
and rjb Bcli'iil — e,jj,,T) transitions. 



2 

max 



and decay widths r<>(in 10"^^ GeV) for Be (D, r/c, B, Bs)£ue 



Mode 



LinearpO] EFG [9, 10] IKS [5] NW [15] HNV [16] AKN [48] CD [14] WSL [17] 



B^^D /+(o^(0) 

/+(l7max) 
/o('7max) 


0.086[0.079] 


0.14 


0.69 


0.1446 


0.089 


0.16 


1.129[0.789] 


1.20 


2.20 


1.017 


0.59 


1.10 


0.673[0.554] 


0.64 








0.59 




0.021[0.014] 


0.019 


0.26 


0.020 


0.005(0.03) 


0.043 




0.019[0.012] 













-Be ^ 77c /+(0)(0) 
/+('i'max) 
/o(gmax) 


0.482[0.546] 


0.47 


0.76 


0.5359 


0.49 


0.622 




0.61 


1.084[1.035] 


1.07 


1.07 


1.034 


1.00 




0.94 


1.10 


0.876[0.872] 


0.92 






0.91 






0.86 




6.93(7.95] 


5.9 


14.0 


6.8 


6.95 


8.6 


2.1(6.9) 


9.81 




2.31[2.46] 




3.52 




2.46 


3.3 ±0.9 







Bc^B /+(o^(0) 

/+(9max) 
/o(q'max) 


0.464[0.428] 


0.39 


0.58 


0.4504 


0.39 


0.362 




0.63 


0.729[0.647] 


0.96 


0.96 


0.6816 


0.70 




0.66 


0.97 


0.572[0.570] 


0.80 






0.71 






0.81 




0.84(0.69] 


0.6 


2.1 


0.638 


0.65 




0.9(1.0) 


1.63 


Tm 


0.80(0.67] 








0.63 









Be ^ Bs /+(o)(0) 

/+('i'max) 
/o('?max) 


0.570(0.574] 


0.50 


0.61 


0.5917 


0.58 


0.564 


0.73 


0.802(0.771] 


0.99 


0.92 


0.8075 


0.86 


0.66 


1.03 


0.685(0.716] 


0.86 






0.86 




0.87 




15.45(15.20] 


12 


29 


12.35 


15.1 


15 11.1(12.9) 


23.45 


Tm 


14.61(14.40] 








14.5 







rib 



Be 



/+(o^(0) 

/+('i'max) 

^ ) 

max/ 



/o(g, 



0.341(0.523] 
0.976(0.918] 
0.811(0.839] 
4.64(7.94] 
1.57(2.11] 



tial model predicts T{B; B^-i) = 22.4 [19.9] eV - 
1836 [1631] eV for Am = 50 MeV - 220 MeV. This 
sensitivity for the B* radiative decay may help in deter- 
mining the mass of B* experimentally and pinning down 
the best phenomenological model. Other magnetic dipole 
decays V — > P7 of various heavy-flavored mesons such as 
{D,D*,Ds,Dl,Tj„J/i,) andjB,B*,B„B*,rjb,T) using 
our LFQM can be found in [28| . 

VII. SUMMARY AND DISCUSSION 

In this work, we investigated the exclusive semileptonic 
Be {D,r]c,B,Bs)£ve, rib B^lvt {I = e,fj.,T) de- 
cays and the magnetic dipole B* Bcj decay using our 
LFQM constrained by the variational principle for the 
QCD motivated effective Hamiltonian with the linear (or 
HO) plus Coulomb interaction. Especially, we obtained 
the new variational parameter f3cb for the bottom-charm 
sector and predicted the mass eigenvalues of the low- 
lying Be and B* states. Our new predictions of Mb^ — 
6459 [6351] MeV obtained from the linear [HO] potential 
model is in agreement with the data, = (6276 ± 4) 

MeV [40], within 3% error. We also predicted the un- 
measured mass of B* as Mb. = 6494 [6496] MeV for the 
linear [HO] potential model. Our model parameters ob- 
tained from the variational principle uniquely determine 
the physical quantities related to the above processes. 
This approach can establish the broader applicability of 



our LFQM to the wider range of hadronic phenomena. 
For instance, our LFQM has been tested extensively in 
the spacelike processes [U, [s^l as well as in the timelike 
exclusive processes such as semileptonic [H, [H, Isl] and 
rare [13] decays of pseudoscalar mesons and the magnetic 
dipole V P-f* decays [Ulli]. 

The weak form factor f±(<f') for the semileptonic de- 
cays between two pseudoscalar mesons and the decay 
form factor Fb^bA^^) for the B* Bcj decay are ob- 
tained in the g+ = frame {q^ = < 0) f-iid then 
analytically continued to the timelike region by chang- 
ing to —q^ in the form factor. The covariance (i.e., 
frame independence) of our model has been checked by 
performing the LF calculation in the 17+ = frame in 
parallel with the manifestly covariant calculation using 
the exactly solvable covariant fermion field theory model 
in (3 -|- l)-dimensions. We found the zero- mode contri- 
bution to the form factor f-{q^) and identified the zero- 
mode operator that is convoluted with the initial and 
final state LF wave functions. We calculated the decay 
constants of (Sc, -B*) mesons and the decay rates for the 
exclusive B^ {D,ric, B, Bs)£i'i and rji, Bcii^e decays 
and compared with other theoretical approaches. Par- 
ticularly, the decay constants for {Bc,B*) mesons and 
the decay rate for rjb Be process are quite sensitive to 
the choice of potential within our LFQM. From the fu- 
ture experimental data on these sensitive processes, one 
may obtain more realistic information on the potential 
between quark and antiquark in the heavy meson sys- 
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tern. 

For the radiative B* B^'y decay, we find that the de- 
cay width r{B* Bcj) is very sensitive to the value of 
Am = AIb' — Mb^- This sensitivity for the B* radiative 
decay may help in determining the mass of B* experimen- 
tally. Since the form factor FwBS'f') ^'^^ the radiative 
B* — > i3c7 decay presented in this work is analogous to 
the vector current form factor g{q^) in the weak decay of 
ground state vector meson to ground state pseudoscalar 
meson, the ability of our model in describing the radia- 
tive decay would therefore be relevant to the applicability 
of our model also for the weak decay. Consideration on 



such exclusive weak decays in our LFQM is underway. 
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